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ABSTRACT 

We construct T 2 -symmetric charged AdS black holes of conformal gravity. The most 
general solution up to an overall conformal factor contains three non-trivial parameters: the 
mass, electric charge and a quantity that can be identified as the massive spin-2 hair. We 
study the Dirac equation for the charged massless spinor in this background. The equation 
can be solved in terms of the general Heun's function for generic frequency u and wave 
number k. This allows us to obtain the analytic Green's function G(oj, k) for both extremal 
and non-extremal black holes. For some special choice of back hole parameters, we find 
that the Green's function reduces to simpler hypergeometric or confluent hypergeometric 
functions. We study the Fermi surfaces associated with the poles of the Green's function 
with vanishing uj. We find examples where the Fermi surfaces for non- Fermi liquids as well 
as the characteristic Fermi ones can arise. We illustrate the non-trivial differences in the 
Green's function and Fermi surfaces between the extremal and non-extremal black holes. 
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1 Introduction 



The AdS/CFT correspondence [HE1E] provides a new method of studying strongly coupled 
field theories by constructing the corresponding gravitational duals which are classical and 
relatively simpler. In particular, the Green's functions in the dual field theories can be 
derived by studying the wave equations of the corresponding bulk fields such as the scalar 
or spinor fields. One of the most studied class of gravitational backgrounds is the charged 
black holes that are asymptotic to anti-de Sitter spacetimes (AdS) whose bulk vectors are 
dual to the currents in the boundary field theories. Such bulk geometry allows one to 
investigate properties of some strongly coupled fermionic system at finite charge density 
such as non- Fermi liquids [4j. Analogous to those of the scalar operators, the fermionic 
Green's function can be obtained by analysing and solving the bulk Dirac equation. 

The procedure of obtaining Fermi surfaces from the AdS/CFT correspondence was 
spelled out in detail in [21 E]. (See also the review [7j-) In general, for an electrically-charged 
black hole background, one can write down the Dirac equation for a charged spinor. The 
mass of the spinor is typically chosen to be zero in practice except for analyzing the qual- 
itative properties. If one knows how to solve the equation, one can read off the Green's 
function in the momentum space from the asymptotic behaviour of the wave solution after 
imposing some appropriate horizon boundary condition. The Green's function is expressed 
in terms of the black hole parameters and the electric charge q of the spinor, as well as 
oj, the frequency and k, the wave number. The poles of the Green's function imply the 
existence of quasi particles. A Fermi surface kp is defined as a pole of the Green's function 
in the momentum space with vanishing u. 

Since the Fermi surfaces are determined by the uj = Green's function, it is not necessary 
to know the Green's function for general (lj, k). It appears sufficient to obtain wave solutions 
with co = and the associated Green's function G(u = 0, k)^ Nevertheless, it is of great 
interest to study the behaviour of the Green's function for small u) on or near the Fermi 
surface. For extremal solution with zero temperature, a procedure was developed in |6] 
to calculate such a Green's function. It takes the form 



where hi,vp,h,2 and v^ F are constants that can be determined. In particular, vf here is 

the Fermi velocity. It was known in the condensed matter physics that a Landau-Fermi 

1 As we shall see later, the situation is subtler for non-extremal backgrounds, where there is no criterium 
to impose the horizon bounary condition on the oj — solutions alone. 




hi 
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surface of Fermi liquids is given by vp = 1. The Green's function with v > ^ shares the 
same characteristics of that near the Landau-Fermi surface. The situation with up < \ is 
very different and the Green's function is associated with some non- Fermi liquids. 

In practice, the above description is optimistic. Typically the Dirac equations cannot be 
solved analytically. In [6] , the extremal Reissner-Nordstrdm AdS black holes in general 
dimensions were used as gravitational backgrounds. It turns out that the Dirac equation 
cannot be solved even for ui = 0. The above procedure were performed numerically. Never- 
theless, it shows that Fermi surfaces can arise from the extremal Reissner-Nordstrdm AdS 
black holes, and Green's functions for small oj near the Fermi surfaces indeed take the form 
([I]) . Furthermore it was shown in [5j [6] that Fermi surfaces for non- Fermi liquids can also 
arise, which are characterised by vp < |. 

Although a physical system can be adequately analysed numerically, it is more satisfying 
if one has analytic results. Reissner-Nordstrdm AdS black holes in D = 5 and D = 4 
are special cases of charged AdS black holes in gauged supergravities. The more general 
solutions are the three-charged [Hj and four-charged black holes [2] respectively, and their 
type IIB and M-theory embedding were given in |lUj . In [TT], the five-dimensional black 
hole were considered and reduced by setting two charges equal while turning off the third 
charge. Taking a certain suitable extremal limit, Fermi surfaces of the resulting G(0, k) 
can be read off analytically. This enables one to construct an analytical expression for 
the constants hi,VF,h2 and Vf, F in ([!]). However, the extremal black hole suffers from 
a curvature singularity that is on the horizon. The general G(u, k) is unknown for this 
example. Moreover, such a Green's function associated with a general non-extremal charged 
black hole is hitherto unknown. On the other hand, the examples considered in [5j [6] and 

are already among the simplest black holes in the usual two-derivative theories of gravity 
or supergravity. It is unlikely to find new black holes in these theories without exotic matter 
so that the Dirac equation becomes exactly solvable. 

Charged black holes can also arise in conformal gravity. In four dimensions, the Weyl- 
squared gravity is conformally invariant. Its minimum coupling to the Maxwell field pre- 
serves the conformal symmetry. Conformal gravity has attracted renewed interest recently. 
It was shown that the Einstein- Weyl gravity has a critical point [12] for which the seemingly 
inevitable ghost massive spin-2 excitations are replaced by modes that fall off logarithmi- 
cally. It was proposed that Einstein gravity can emerge from conformal gravity in the 
infrared region |13| . Furthermore, conformal gravity can be supersymmetrised in the off- 
shell formalism |14| . 



4 



The allowance of a vector field in conformal gravity implies that the theory should 
have a variety application in the AdS/CFT correspondence. It can provide gravitational 
dual description of certain strongly interacting fermionic system at finite charge density, on 
which we focus in this paper. The most general spherically-symmetric black holes carrying 
electric charges, up to an overall conformal factor, were obtained in |15j . For our purpose 
in this paper, we construct the analogous solution with the torus (T 2 ) topology. We study 
the charged Dirac equation in this background and find that it can be solved exactly for 
general (u,k). This enable us to construct the general Green's function G(u,k) for non- 
extremal black holes as well as their extremal limits. The results are expressed in terms of 
general Heun's functions, whose properties are less known. We thus consider some special 
sub-classes of the black holes for which the wave solutions are reduced to hypergeometric 
functions. This enables us to study the Green's functions in greater details, not only for 
extremal black holes but also for non-extremal ones. 

The paper is organised as follows. In section 2, we review conformal gravity and the 
spherically-symmetric black holes. We then construct the static black holes with torus T 2 
and hyperbolic H 2 horizons. The most general solutions contain three non-trivial parame- 
ters, in addition to the cosmological constant. We study the global structure and analyse the 
thermodynamics. We then consider two special subclasses of solutions with two parameters. 
One was obtained previously in [16J and it can be viewed as the pseudo-supersymmetric 
solution in the corresponding conformal off-shell supergravity. The two-parameter solution 
has both inner and outer horizons. The solution is in general non-extremal, and it becomes 
extremal when the two horizons coalesce. The other class is the general extremal solution 
which has also two non-trivial parameters. In section 3, we consider Dirac equations for 
a charged massless spinor in the black hole backgrounds. For the "BPS" black holes, the 
wave equation can be solved exactly. We thus obtain the analytic Green's function for 
general (ui, k); it can be expressed in terms of hypergeometric functions. We also obtain the 
Green's function in the extremal limit, and it turns out be expressible in terms of confluent 
hypergeometric functions. We study the Green's functions G(0, k) and obtain the Fermi 
surfaces. The results allow us to study the intrinsic differences on the boundary field theory 
between the extremal and non-extremal backgrounds. We accomplish the above in section 
4. 

In section 5, we consider the two-parameter family of extremal black holes, and we 
find that G(0, k) can be expressed in terms of hypergeometric functions. This allows us to 
determine the Fermi surfaces. In section 6, we demonstrate that the Dirac equation can 
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be solved even in the most general black hole backgrounds. However, the resulting Green's 
function G(u, k) is expressed in terms of the general Heun's functions. We discuss how the 
result can be reduced to the previous simpler examples. Owing to the complication of the 
Heun's function, it is difficult to discuss many subtle properties of the Green's function; 
nevertheless, we are able to find many examples of Fermi surfaces. The Green's function 
has a fascinating rich structure of spiked maxima in the (oj, k) plane. We conclude our 
paper in section 7. 

2 Charged black hole in conformal gravity 
2.1 Review of conformal gravity 

Analytical charged AdS black holes also arise in conformal gravity in four dimensions. 
Conformal pure gravity is constructed from the Weyl-squared term. Its conformal symmetry 
is preserved when it couples minimally to the Maxwell field. The Lagrangian is given by 

er l C = \aC» upa C^ pa + \aF 2 , (2) 

where e = y/—g, F = dA and C pupa is the Weyl tensor. Note that the absolute value of the 
coupling of the Maxwell field can be arbitrary defined by a constant scaling of A. However its 
sign choice is non-trivial. We have made this selection with the following two considerations. 
The first is inspired by critical gravity |12j . If we add an cosmological Einstein Hilbert term 
e (-R + 6) to the Lagrangian, the critical point is precisely a = — in which case, the vector 
A is non-ghost-like for our choice. The second is the consideration that Einstein gravity 
may emerge from conformal gravity in the infrared limit [13] . The conformal gravity with 
a = — j can give rise to Einstein gravity e(R + 6). Having determined the theory, it is 
straightforward to derive the equations of motion, given by 

V^ = 0, -a(2W CT + R^Cppvv - la(F% - \F 2 g pv ) = . (3) 

It is worth pointing out that conformal gravity can be supresymmetrized in the off- 
shell formalism. In J\f = 1, D = 4 off-shell supergravity, the bosonic fields consist of 
the metric, a vector and a complex scalar S + iP. The fermionic field involves only the 
off-shell gravitino W Up to and including the quadratic order in curvature, the theory 
allows four independent super-invariants. These comprise a "cosmological term |17j." the 
Einstein-Hilbert term \1<$\ 119] , and two quadratic-curvature terms [T3], one formed using 
the square of the Weyl tensor, and the other formed using the square of the Ricci scalar. 
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Supersymmetric solutions in these higher-derivative off-shell supergravities were obtained 
including Lifshitz, Schrodinger, gyratons [TB I |20~ 1 [2T | [2"2"j. 

The Lagrangian ^ turns out to be the bosonic part of the super invariant, except 
that an analytical continuation of A — > iA was performed. The resulting theory is called 
pseudo-supergravity with the Killing spinor equation [TBI [22] : 

5% = -D^e + l(2A„ - r^A u )T 5 e - |r M (S + r 5 P)e = . (4) 
2.2 General black holes and thermodynamics 

Up to an overall conformal factor, we find that conformal gravity admits the following static 
charged black hole solution^] 

ds 2 = _ fdt 2 + ^l +r 2 dn 2 A =-®dt, 

f ' r 

f = -|Ar 2 + cir + c + - , 3ad + e 2 + Q 2 = c\ , (5) 

where e = 1,0, —1 for df2| e being the metric for the unit sphere S 2 , torus T 2 and hyper- 
bolic H 2 . The solution for e = 1 was given in |15j . Note that although ^ contains the 
Schwarzschild black hole as a special solution, it does not contain the Reissner-Nordstrom 
black hole. This is important for this paper since the Dirac equation in the Reissner- 
Nordstrom black hole is not analytically solvable. 

We now study the global structure and thermodynamics of the general black hole. (The 
thermodynamical properties for Q = were given in [23].) Let r$ > denote the largest real 
root of /(r). The metric on and outside the horizon runs from r = r$ to the asymptotic AdS 
boundary at r = oo. It is advantageous to express Q in terms of the remaining parameters 
and ro, given by 

Q = Jc 2 -s* + 3 ^- + 3 ^-Ar d. (6) 

The temperature can be obtained by the standard method of Euclideanising the metric and 
demanding the regularity on the horizon. The entropy can be calculated directly from the 
Ward formula. They are given by 

6d + 3c r + Ar[j 1 a(3d + c ro)w 2 M , 

T = — g > S = 6 aea; 2 ~ R • (7) 

127T7q 6r 



2 As we shall see later, the wave equation is invariant under the conformal transformation that is r- 
dependent only. It follows that the Green's function and Fermi surfaces are independent of the T 2 -symmetric 
conformal transformation. 
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The electric potential and charge are 

* = Qe=^. (8) 

r 12tt V 7 

As discussed in [23], it is also necessary to introduce a pair of thermodynamical quantities 
that can be viewed as the massive spin-2 hair: 

The free-energy can be obtained from the Euclidean action, given by 



, , _, /•oo , . ai 2(cn — e)ern + (3e - 

F - £ ^(ii^P + H - 1 2 J ^ • ^ 

Note that for the T 2 topology e = 0, the sign of the free energy is determined by the sign 
of the parameter d. We can also treat the cosmological constant A as a thermodynamical 
variable, since it is an integration constant in the solution. Its thermodynamical conjugate 
is given by 

° = w- <»> 

The energy calculation is more sutble. It was shown [23] that both the ADT |24[ [25] and 
AMD |26[ [2?1 [28"1 [29] methods would give rise to divergent results for this type of black 
holes. It can however be calculated using the method outlined in [23] by employing the 
Noether charge associated with the time- like Killing vector. We find 

a(c - e)(Arg - 3c ) a(2Arg - c + e)d 
E= 72^o + 24^g • ( 12 ) 

It is then straightforward to verify that the following thermodynamical relations are satis- 
fied: 

dE = TdS + <S>dQ e + ^dH + 6dA , F = E — TS — $Q e . (13) 
It is worth pointing out that the Smarr formula now becomes 

E = 2GA + ^E. (14) 
2.3 A class of "BPS" black holes 

The general solution ^ has a total of four parameters. In this paper, we shall concentrate 
only on the solutions that are asymptotic to the AdS. It follows that the parameter A is 
trivial and can always be set to a fixed negative constant. We shall let it be A = —3 so that 
the AdS has the unit radius. The general solution has then three non-trivial parameters, 
associated with the mass, charge and the massive spin-2 hair. 
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As discussed in section 2.1, conformal gravity can be supersymmetrised in the N = 1 off- 
shell formalism. Although in our Lagrangian ^ the vector field has changed from A to \A 
compared to the supergravity theory, the resulting Killing spinor equation Q is still valid. 
It was shown in |16j that for the static black hole configuration, the existence of a Killing 
spinor requires that e = 0, i.e. the topology of the black hole horizon is T 2 . Furthermore, 
we must have d = in Thus the "BPS" solutions take the form [16] : 

ds 2 = -fdt 2 + — + r 2 {dx 2 + dy 2 ), A = adt , (15) 

where /, (f> are functions of r only: 

f = r 2 -fir + Q, a = -®, (16) 

r 

The solution has a local Killing spinor of Q provided that 

S= ^~^ + 6r 2 

2rVf K ' 

It is convenient to parameterise the solution as follows 

/ = (r — r_)(r — r+) , Q = r_r + , = r_ + r + , r + > , r + > r_ . (18) 

Thus, in general the solution describes an asymptotic AdS charged black hole with the outer 
horizon at r = r + ; r = r_ is the inner horizon. The thermodynamical quantities are 

r i — r_ 



S = —QOtr-r + uj2 , 



ar_r + w 2 

^ = r_+r + , *= — , 

ar-r + (r +r + )a^ 

24tt v 7 

They satisfy the following thermodynamical relationship 

dE = TdS + $dQ e + VdE , F = E-TS- $Q e . (20) 
For this subclass of solutions, the free energy F vanishes identically. As we may see from 



the expression of the free energy (10) that for e = 0, the sign choice of the parameter d 
dictates the sign of F. Thus the above "BPS" black hole sits in the critical point of the 
Hawking-Page global phase transition. The internal energy at this phase transition point 
is E = TS + $Q e - 

It is worth pointing out that although there exists a local Killing spinor for general 
parameters, the function S blows up on the horizon r = r+, implying that the solution 
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is not truly supersymmetric. However, this defect disappears in the extremal limit with 
r + = r_ = 7~o, giving 

ds 2 = _ (r _ ro) 2 dt 2 + dr 2 + ^ 2 + 2 A = r (r-r ) dt 

(r — r^y r 

In this limit, we have T = = F and E = &Q e . Note that if we require that S be constant, 

there exists an alternative extremal solution: 

ds 2 = - ^^dt 2 + + r 2 (d X 2 + dy 2 ) , A = ^ ' 7 ^ dt . (22) 

Both extremal solutions run from the AdS2 x T 2 horizon at r = to AdS4 at r = oo. 
It was shown that these two extremal solutions are related one to another by an overall 
conformal scaling [16J. A priori, in this paper they should be treated different solutions 
since the coupling of conformal gravity to a free scalar or spinor breaks the conformal 
symmetry. However, as we shall see later, Fermi surfaces are invariant under the conformal 
transformations that preserve the T 2 isometry. 

2.4 Classes of non-BPS black holes 

The general solution ([5]), which has three non-trivial parameters, is rather complicated for 
our purpose of studying the dual conformal field theory in some detail. The previous "BPS" 
limit is a way of simplifying the solution, corresponding to setting d = 0. We consider three 
more simplified cases with non-vanishing charge. The first corresponds to c\ = 0, the second 
corresponds to cq = and the third is the general extremal solution. In all these three cases, 
there are two non-trivial parameters, one of which is the electric charge. For our purpose 
later, we shall consider only the T 2 topology, namely s = 0. Without loss of generality, we 
set A = -3. 

Case 1: c\ = 

The solution is given by 

f = r 2 + Q + - r , a = -£. (23) 

It admits an extremal limit, for which / has a double zero at the horizon r = ro- The 
parameters and the solution are given by 

Q = -3r 2 , d=2rl f = {r - r ) 2 {r + 2r») ^ g = _ 3r (r - r ) 

r r 

Note that we have made a gauge choice so that a vanishes at r = ro. It is of interest to 

note that the same local metric also describe a non-extremal black hole when r is negative. 
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(The horizon is at r = — 2tq and the asymptotic AdS arises as r — > — oo.) Alternatively, we 
can take a view that non-extremal black holes arise when tq becomes negative. Thus the 
solution describes a black hole for any real and non- vanishing tq. 

Case 2: Co = 

The solution is given by 

f = r 2 + Cl r-^, a = -®. (25) 
6c\r r 

It also admits an extremal limit. The corresponding parameters and the solution are given 
by 



3 , (r - r ) 2 (2r + r ) 3r (r - r ) 

* z 2r 2r 



3^2 



As usual, we have made a gauge choice so that a vanishes at r = r$. As in the case 1, the 
same local metric also describes a non-extremal black hole when r is negative. (The horizon 
is at r = —tq/2 and the asymptotic AdS arises as r — > — oo.) 

Case 3: General extremal solution 

It follows from ([5]) that the function / for the general extremal solution must take the 
form f = (r — ro) 2 (r — r\)/r with ro > and ro > r%. This implies that 

co = r (r + 2n) , ci = -2r - n , d = -r§n , Q = r (r - n) . (27) 

This leads to the general extremal solution 

(r - r ) 2 (r - n) {r - r )(r - n) 
J — > « — • l^oj 

(Note that if we let ri > ro, the solution becomes non-extremal black hole with horizon 
located at r = r%.) If we set r\ = 0, the solution becomes the "BPS" extremal solution 



(21). If instead we let r\ = 2tq or r\ = ^ro, we obtain the extremal solutions (24) and (26) 
respectively. 



3 Charged spinor and the Green's function 

In the previous section, we have constructed a variety of charged black holes in conformal 
gravity that are asymptotically AdS. In this section, we review and discuss the formalism 
of deriving the Green's function from a charged spionor satisfying the Dirac equation [51 E]. 
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The Dirac equation for a charged massless spinor tp in a generic gravitational background 
is given by 

/(^^t-i^^O. (29) 

Massive charged Dirac fermions were also considered in 0(6]; however, they tend to lead to 
unsolvable wave equations and are necessary mainly for discussing some qualitative features. 
We shall not consider them in this paper. To derive the formalism, we consider a more 
general background in the form of 

dr 2 

ds 2 = -fdt 2 + — +r 2 {dx 2 + dy 2 ) , A = adt. (30) 

Note that the vector has a gauge symmetry A — > A + dA. In the Dirac equation, this gauge 
symmetry has the effect of creating a phase factor in the fermion field ip. In particular, the 
constant shift in At implies a constant shift in uj. Thus in this paper, we shall always make 
a gauge choice such that A vanishes on the horizon by an appropriate constant shift in At. 
This gauge choice makes the boundary condition of the wave solution on the horizon easier 
to define. 



The vielbein and the corresponding spin connection for ( 30 ) are given by 

e° = /2 dt , e 1 = rdx , e 2 = rdy , e 3 = —j= , 

yh 

u° 3 = \h?f'f- x e° , w\ = h^r~ l e l , u 2 3 = hkr^e 2 . (31) 



The contribution of the spin connection to the Dirac equation (29) can be absorbed by the 
scaling of the field \E' = (—gg rr )^^> . In general a wave function involves the frequency to 
and the wave numbers k x and k y . Owing to the symmetry in the (x,y) plane, we can set 
k y = without loss of generality. In the Fourier mode \P ~ e -^t+ikx^, w j^] 1 ^ e momentum 
lying in the x direction only, the Dirac equation becomes 



- if- a (u + #) 7 U + hi^dr + ir^kj 1 ) * = . (32) 
Adopting the following conventions for the gamma matrices 



( itfi 


0^ 


i M 


1° 


io-i j 








7 3 = r 7 s 



icr 2 
-icr 2 



the Dirac equation reduces to two decoupled equations 

+ q<f>)*! + h^a 3 d r + (-l) a r -1 iA:<72]&« = , (33) 
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where a = 1,2. Note that only three gamma matrices (7°,7 3 ,7 ) appear in the wave 



equation ( 32 ) . The detail gamma matrix decomposition above is unnecessary except for the 
three. Thus, it is rather straightforward to generalise the formalism to arbitrary dimensions 

61- 

The equation for ifj2 is related to the equation for ipi by k — > —k. Each ip a is a two- 
component spinor. Let ip\ = {u\,U2) T and u± = u\ ± iu2, we find that 

u' + + Xi(r)u + = A 2 (r)«_ , u'_ + Ai(r)n_ = \ 2 {r)u + , (34) 

where 

i(w + ga) i/c 

Al ( r ) = 7f^' A 2( r ) = 7=- (35) 



(Note that if we had considered the massive Dirac equation, we would have instead A 2 (r) = 
VK — ^7^0 Here the bar denotes the complex conjugation. It follows that 

u '+ + Pi ( r ) u + + = , (36) 

u"_ + Pl (r)«'_ + P2(r)«_ = , (37) 

where 

Pl(r) = -^, P2(r) = |Ai| 2 -|A 2 | 2 +pi(r)Ai + A , 1 . (38) 

The advantage of writing the decoupled equations for u± instead of Ui is that these equations 
now involve only the metric functions rather than their square roots. In this paper, we opt 



to solve for u_ first from the second-order differential equation (37). We then obtain u + 



from (34), namely 

v! + \\U- , 
u+ = ~ . (39) 



Note that although from (36) u + simply takes the form that is complex conjugate of U-, the 



integration constants can only be determined through (39), so that the full set of solutions 



for (u+,U-) have two independent integration constants rather than four. 



As we shall see in the next section, for the "BPS" solutions (18), the wave functions u+ 
and U- can be solved analytically for general u and k, in terms of hypergeometric functions. 
(The analysis on the general black hole is given in section 6.) By choosing the appropriate 
boundary condition on the horizon such that the waves are in-falling only, we can then read 
off the Green's function by studying the behaviour of the wave functions at AdS boundary, 
namely 

G 1 = -Go 1 =lim^ = -i lim U+ - — . (40) 

r->oo u\ r->oo U+ + U_ 

13 



Here, G\ and G2 are diagonal entries of the 2 by 2 matrix of the Green's function associated 
with the massless spinor [5l|6l[TT]. In particular G\ corresponds to the right-handed spinorial 
operator and the G2 corresponds to the left, which we shall not consider in this paper. In 
this "standard" quantization, the function u\ is treated as source and 112 is the response. 
It is clear that u\ and 112 are symmetric from the point of view of the wave function. This 
leads to an alternative quantization in which the roles of u\ and 112 are reversed. In this 
case, the Green's function of the right-handed spinor is given by — G2 instead. 



For a generic black hole of the form ( 30 ) , if we could obtain the general wave solution 
and hence the Green's function, we could then determine the Fermi surfaces k = &p which 
are defined as poles of the Green's function of vanishing oj. In practice, However, the Dirac 
equation cannot be solved exactly for generic (ui,k). One can then try to solve the Dirac 
equation for simpler case with u = 0. This is enough for obtaining G(0, k) and determining 
a Fermi surface. It is still of interest to study the behaviour of the Green's function for small 
uj on or near the Fermi surface. In order to derive such properties on the Fermi surface, 
a procedure was developed in [SJ E] for the extremal Reissner-Nordstrom black holes in 
general dimensions. Unfortunately, there was no analytic solution even for uj = 0. The 
above procedure was only proceeded numerically. Although the the Green's function ([!]) 
can be reproduced numerically even for non- Fermi liquids, it is more satisfying to obtain 
some analytic expression for such a Green's function. 

In |llj . an example of extremal background was discussed in which G(0, k) and hence 
the Fermi surfaces could be obtained exactly. This gives rise to an analytical expression 
for the Green's function of small uj, namely Q. However, the general expression of G(uj, k) 
is still unknown. In addition, an analytic expression of Green's function for non-extremal 
black hole is also lacking in the literature. This story changes if we use instead the black 
hole backgrounds of conformal gravity: the Dirac equations become exactly solvable! We 
shall proceed in the next section. 

To end this section, we would like to comment on the effect of the conformal symmetry 
on the Fermi surfaces. If a conformal transformation does not alter the global structure 
such as the asymptotic infinity and the horizon, the resulting new solution can be viewed 
as equivalent to the original one. However, the Dirac equation is not in general invariant 
under the conformal transformation. If we restrict the conformal transformations to be r- 



dependent only, the equations in (34) for u± are invariant. Thus the Green's function G(uj, k) 
and hence the Fermi surfaces are invariant under this class of conformal transformations. 



14 



4 Exact Green's function from the "BPS" black holes 



As we shall demonstrate in section 6, the Dirac equation (29) on the general T 2 -symmetric 



background ^ with e = can be solved exactly, and the analytic Green's function for 
general (oo, k) can be obtained. However, the results are expressed in terms of the general 
Heun's functions, which are less studied. In fact there is no Heun's function defined in 
Mathematica, and very limited properties of these functions are coded in Maple. In this 
section, we study the Dirac equation in the simpler "BPS" backgrounds, given in section 



2.2, namely, the black holes given by (15) with (18). We make a gauge choice 



A = (r_- T ^)dt, 



(41) 



so that it vanishes on the horizon r = r + . 



4.1 Green's function from Non-extremal black holes 



Let us first consider the general solution (18) with two parameters (r + ,r_). We have 
.(u) + qrJ)r — qr + r- . , , ik 



Ai(r) = i- 



A 2 (r) 



(42) 



r(r ~ r+)(r - r_) ' ~ r ^J(r - r+)(r - r_) ' 

It turns out that the wave equations for u + and u_ are exactly solvable. We organize the 
solutions as in-falling and out-going modes on the horizon. The in-falling solutions are given 

by 



2k 



Aoo + i(r+ — r 
2F1 



1 — iq + v, 1 — iq — u; | 



2iO; 



r_(r — r_) 



11 _ 



r iq (r — r + ) lfi (r — r_ 



-iq+iQ , 



r(r + — r_) 



where 







2-^1 



r+ — r_ 



ig + 1/, — iq — u; \ — 2if2; 



T—\r — Ta 



r(r+ — r_) 



CO 



Ik 2 
r+r- 



Q 



(43) 



(44) 



Note that T is the temperature of the background black hole, given in (19). These solutions 



are called in-falling since near the horizon, the wave function behaves like 



tp ~ exp f — iwt + iil log(- 



4vrT 



The outgoing wave functions are given by 



out 



U 



out 



(45) 



(46) 
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where * denotes the complex conjugate. The outgoing solution should be discarded for the 
black hole background. Of course, if we reverse the time, the outgoing solution becomes 
in-faling. Thus the final answer of the Green's function is independent of whether we choose 
in- falling or outgoing solutions. 

Following the discussion in section 3, we find that the Green's function for general oj 
and k is given by 



where 



7 



2 .Fi[l - v - iq, 1 + v - iq; § + ifi; -j^E^] k 
2 Fi[-i> -iq,u- iq; \ + ifi; ~^5pz] (r+ - r_] 
2 Fi[l-u-iq, 1 + v - iq; § + iO; k 



2 Fi[-p -iq,v- iq; \ + iU; (4vrT) 



(48) 



The Green's function for to = is then simply given by (47) and (48) with Q set to zero. 
It is clear that the general Green's function G(u, k) is complex and it remains complex for 
ui = 0. By contrast, as we shall see in the next subsection, G(0, k) for extremal black hole 
becomes real. Interestingly, in non-extremal backgrounds, the quantity v does not have any 
particular physical importance and there is no requirement that it be real. The situation is 
changed in the extremal limit, where v has to be real for stability. 

We find that for some special parameter choices, the hypergeometric functions in the 
Green's function degenerate. For example, if we let q = = cj, we have 

i + tanh ( , 2k arcsinh- 



G = i Vyf+E [±ZIW . (49 ) 



tanh ( , 2k arcsinh — ) 



For v = \ and to = 0, we have 



where 



G = i 2 ~ i V / ( 1 + 16 g 2 ) r + r -^ ( 50 ) 
2 + i^(l + 16g 2 )r+r_ 7 



2i V V^T^ 



1 + i 



! = r , — '-T7ZZ — • (51) 



(1 - 4ig) A /r_(r+ - r_) A _ . 



+ 1 + i ' ' ^ " 



r_|_— r_ / V V r + — r 



4.2 Green's function from the extremal black hole 

We now consider the extremal case with r_ = r + = ro- This case was discussed in some 
detail in [30] . Although the results in this case can be obtained by taking a subtle extremal 
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limit on the hypergeometric functions in the previous subsection, it is of interest to study 
the extremal case on its own and we shall discuss this example independently. In this case, 
we have 

(u + qr )r - qr 2 ife 

Ai = 1 , \ 2 = -— - (52J 

r(r — r§y r(r — r$) 

The in-falling solutions of u± can be obtained, given by 

u _(r) = z~*W K , v (z) , u + (r) = -z^W-^z) , (53) 

where W K ^(z) denotes the Whittaker function. The new notations above are specified as 
follows 

K=±+iq, v = . / K - q 2 , z= . , 54 

V r {r ~ r o) 

and k* is the complex conjugate of k. As we shall discuss at the beginning of subsection 
4.4, the quantity v has to be real in this extremal limit, which implies that \k\ > \qro\. The 
sign choice of v is however conventional, and it has an effect on defining which mode is 
in-falling. The following identities are useful 

= (k-Iz)W k>v (z)-(^-(k-1) 2 )W k - 1>v (z) 

= -(K-\z)W K , v -W K+lyU {z). (55) 

As r approaches the horizon r = ro, we have z — > oo, and the Whittaker function behaves 
as 

W KtU ^e-^ z z K {l + 0{z- 1 )). (56) 

It is then straightforward to recognise that ip is in-falling only. The retarded Green's 
function, associated with the in-falling modes, can now be read off straightforwardly, given 
by 

Gi = -G~ 2 l = i^ , (57) 
1 + 7 

where 

7 = -V^ 1 • (58) 

r U(-iq + v,l + 2v,- 2 %) 

Here we have used the identity 

W KtV (z) = eM-h)^U{u-k + \,l + 2u-z), (59) 
where U is the confluent hypergeometric function of the second kind, defined by 

U ( a > 6 ' z ) = r(a-6 + l) li?l(Q; b ' Z) + ^T(6T^ 6li?l(a -b+l;2-b;z). (60) 
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We can also obtain the result directly from (47) and (48) by taking a limit of r_ — > 
r+ = tq. Since the last argument in the hypergeometric function in (48) blows in this limit, 
it is convenient to use the following identity 

2 F 1 (a,b;c;z) = - °j {-z)- a 2 F x (a, 1 + a - c; 1 + a - b; z' 1 ) 

1 (ojl (c — a) 

+ l\ C) ^~^ (-z)- b 2 F 1 (b, l + b-c;l + b-a; z" 1 ) , (61) 
1 (a)l (c — b) 

for z (0,1). 

It is worth emphasizing an important difference between the wave functions in the 
extremal and non-extremal black holes. In the extremal background, the general wave 
solution for the u_ with to = is reduced to 

(62, 

(u + is analogous.) This implies that for v > 0, the wave function is well defined on the 
horizon only when c 2 = 0. This allows one to pick the correct mode even when we do not 
have the more general solution with non- vanishing to, in which case, the concept of in-falling 
and outgoing becomes irrelevant. This implies that in the extremal case, the Fermi surfaces 
can be determined without any knowledge with to ^ 0. The situation is very different in the 
non-extremal solution. As we can see in section 4.1, when cj = 0, both solutions of u- are 
well defined on the horizon. We are able to select the right mode only because we have the 
general solution with non- vanishing to, which allows us to separate the in- falling mode from 
the outgoing mode. If we can solve the Dirac equation only for to = in the non-extremal 
background, we shall not be able to fix the boundary condition of the to = solution alone 
on the horizon. Further distinguishing features of extremal and non-extremal backgrounds 
will be addressed in subsection 4.4. 

4.3 Fermi surface from the extremal black hole 

Having obtained the Green's function in the momentum space for general to and k, we can 
now study the Fermi surfaces. They are characterised by the poles of the Green's function 
for vanishing to. Since we have \F\[a,b, 0] = 1 for generic a and b, it follows that for any 



positive v, the second term in the confluent hypergeometric function (60) dominates. Thus 
in the to — > limit, we have 

G(0,t) = l + — £ r = ,/P#, (63) 

utq + i(k - qr ) yk-qr 



18 



which is a real quantity provided that v is real. Note that for half integer u, the identity 



(60) diverges. However, the ratio 7 (58) remains finite and the resulting Green's function 
takes the same form as the above. Thus for a given extremal black hole with vq and the 
given charge q of the spinor, there is only one Fermi surface: k = qr$ > for the standard 
quantization or k = —qr$ < for the alternative quantization, both correspond to the 
vanishing of v. 

Since the Fermi surface in this case occurs at v = only, there is a delicate limiting 
procedure for investigating the Green's function at the vicinity of the Fermi surface. Let us 
first examine the standard quantization case whose Fermi surface occurs at k = kp = qrQ. 
If ijo/rQ approaches zero much faster than (k — kp), we find that 

G X V2kHk ~ kF) . (64) 

Comparing with the general formula of the Green's function near the Fermi surface ([!]), we 
find that 



hi = - y / 2k F (k - k F ) + ■ ■ ■ , h 2 = 0-2(k-k F ) + --- . (65) 

For the alternative quantization, the Fermi surface occurs at k = kp = —qr$ < 0. Again 
in the limit where oj/tq approaches zero much faster than {k — kp), we have 



G 2 = V2k F (k-kp) (66) 

(* " k F ) - qy/ + 2(k - k F )(-^: 
Comparing with the general formula 0, we have 



hi = + v / 2k F (k - k F ) + • • • , h 2 = 0- 2(fc - k F ) H . (67) 

It is instructive also to consider how the Green's function behaves with small 00, after 
we have literally fixed k = kp. We find 

k = qr : Gi = -i + 2q( log( — — ) - 27 - x(i<?) 

V r / 

k = -qr : G 2 = -i + 2g(log(-— ) - 2 7 - X (k)) , (68) 

where 7 is the Euler number and x is the digamma function. In both cases, the divergences 
are logarithmic in small u. This is different from the general formula ([!]) with k = kp. 
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The v = Fermi surface describes an extreme situation of some Non- Fermi liquids. It 
was shown to exist in the extremal Reissner-Nordstrom black hole [5J E] . As we shall see 
later, it can also arise in the more general two-parameter extremal solutions. The example 
discussed in this subsection is somewhat unusual in that this is the only Fermi surface. 



4.4 Imaginary v and Fermi surfaces from non-extremal black holes 

In the above discussion of extremal background, we have considered the case where u, whose 



expression is given by (54), is a real number. The general expression of v for the non- 



extremal solutions is given by (44). Naively, one would expect that this quantity continues 
to be real, in which case, we find that there is no sign of Fermi surfaces in the non-extremal 
background. 

We now consider the possibility that v is pure imaginary. There are two situations that 
this can arise. The first is when r_ is also positive. In this case, v becomes imaginary if we 
have 

k 2 

q 2 >—. (69) 

r + r_ 

The second situation is when r_ is negative. In this case, v is pure imaginary for any (k, q) 
that do not vanish simultaneously. In either case, we shall define v = i/i where 



k 2 

\<?-rlr..„- (70) 
For the extremal black hole, the solution for the wave function u_ with ui = is given 

by 

Thus, we see that if v is real and positive, the wave function is well-defined on the horizon 
by setting C2 = 0. This is in fact exactly the near- horizon behaviour of the general in-falling 



solution (53) in the co — > limit. If on the other hand, v = i/x is pure imaginary, the 
wave function becomes oscillatory on the horizon. The oscillatory mode with u = on 
the horizon makes it impossible to fix the horizon boundary condition. It was argued that 
this corresponds to an infinite number of particle creation and suggests that the geometry 
becomes unstable [6j. Further discussion will be given later. 

This instability does not occur in the non-extremal solutions. It is clear that the hyper- 



geometric functions in ( 43 ) are well behaved when r — > r+. It follows that the fermionic wave 
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function in the near- horizon behaves simply as (45) and there is no oscillatory behaviour 
when oj = 0. Thus as long as the black hole is non-extremal (including the near-extremal 
case), imaginary v does not imply instability and hence should be considered. It is worth 
drawing attention that for an extremal black hole, there usually exists a decoupling limit in 
which the metric becomes its near- horizon geometry AdS2 X T 2 . The quantity v has a clear 
physical interpretation that it measures the conformal weights of the dual operators in the 
boundary field theory of the AdS2- This implies that v must be real. Such a decoupling 
limit is absent in non-extremal solutions and v has no apparent physical significance and 
it could be real or imaginary. In fact, as we see in section 6, in the wave solution and the 
Green's function for the most general black holes, there is no parameter combination that 
can be recognised as v at all. Only in the extremal limit, such a combination appears. The 
appearance of the i/-like quantity in this non-extremal case is accidental. As we have men- 
tioned, there is no Fermi surface at all for real v in non-extremal solutions. Interestingly, 
Fermi surfaces emerge quite frequently when v is taken to be pure imaginary. In other 
words, there is a restriction on the parameter space (k, q) such as \k\ > |qro| in the extremal 
black hole, but such a restriction disappears in the non-extremal background. 
For oj = and v = i/j>, the Green's function is given by 

o<M) = i!=*r, 7= ^•' 1 - i <" + ^ 1 + i ";-^i ; ^T' t . (72) 

Note an important distinguishing feature that the Green's function G(0, k) above is not 
real, whilst it is always real in the extremal limit. This is related to the fact that in the 
non-extremal case, both wave functions are well defined on the horizon when u = 0, and 
the wave functions u+ and U- are organized in terms as in-falling or outgoing modes. This 
implies that ii+ and u_ are not complex conjugate to each other, and the resulting Green's 
function is therefore not real, even when u = 0. 

Fermi surfaces can then emerge when G(0, k) diverges or vanishes (for alternative quan- 
tization) for certain choice of (k, q). Since G(0, k) is complex, it is unlikely that it diverges 
literally, since it would demand the vanishing of the real and complex parts of the inverse of 
the Green's function simultaneously. Our explicit examples however demonstrate that there 
exist values of k for which |G(0, k)\ becomes large in many orders of magnitude. These can 
be regarded as Fermi surfaces in all practical purpose. Since an analytical expression for 
such conditions of (A;, q) is unlikely to exist, we shall present two illustrative examples. In 



both examples we examine the denominator of the Green's function (72), namely 

V = 1 + 2i 7 . (73) 
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The "zeros" of T> signal the Fermi surfaces. Note that in general T> is a complex number. 
For a given black hole with r± fixed, and given q, we have only one real variable k to 
vary. It is hence rather unlikely that both the real and imaginary parts of T> approach zero 
simultaneously for certain values of k. We are looking for surfaces for which the minimum 
of 1 1 + 2vj\ is smaller in many orders of magnitude. 

We first consider an example with both r± being positive. To be specific, let us choose 
r + = 2 and r_ = \. Thus the we have /x = y / q 2 — k 2 . For a given q, k runs from to q. 
For small q, we find that there is no Fermi surface at all. The left plot of Fig. [T] shows both 
the real and imaginary parts of T>, as a function of k, for q = 5. It is clear that although 
the imaginary part of T> passes through a zero, the real part is never particularly small. 
The situation changes for larger q. The right plot of Fig. [T] has q = 10, and an approximate 
Fermi surface emerges at k = 3.633657. 




Figure 1: These are plots of real (blue line) and imaginary (red line) part of T> as functions 
of k for some fixed q. In the left figure, we have chosen q = 5; in the right figure, we have 
q = 10 and a Fermi surface emerges at k = 3.633657. The black hole parameters are r + = 2 
and r_ = 1/2. 

Let us increase the value of q, and we find that more and more Fermi surfaces emerge. 
In Fig. [2j we present the plots for q = 30 and 50. 

It is worth remarking in the above graphs, when the imaginary (red) lines cross zeros, 
the real (blue) lines some times approaches a local minimum which is very small, but not 
necessarily zero literally. For example, for q = 10, the Fermi surface occurs at k ~ 3.633657 
and we have Re(£>) ~ 0.0058. For q = 50, at fei, we have Re(£>) ~ 0.0022; at k 2 , we 
have Re(P) ~ 2.5 x 10~ 6 ; at k%, we have Re("D) ~ 3 x 10~ 10 ; around at k&, we have 
T> ~ —2.2 x 10~ 16 + 1.9 x 10 13 i. Note that in this case, both real and imaginary parts cross 
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Figure 2: These are plots of real (blue line) and imaginary (red line) part of T> as functions of 
k for some fixed q. In the left plot, we have chosen q = 30; in the right plot, we have q = 50 
and four Fermi surfaces occur at fci = 9.5042069087, k 2 = 13.8418437, k 3 = 17.998011219, 
and &4 = 22.945804499291. The black hole parameters are r+ = 2 and r_ = 1/2. 

zero around k^. Thus owing to the fact that T> is complex, the \D\ is not going to hit zero 
literally; however, the local minimums of \D\ can be so small that these points should be 
viewed as Fermi surfaces. In Fig. K3l we plot the absolute value \D\ for the q = 50 example. 




10 20 30 40 



Figure 3: This is a plot of \T>\ with respect to k, for q = 50. The black hole parameters are 
r + = 2 and r_ = 1/2. 

We have presented the effective Fermi surfaces so far for the given q's. It is instructive 
to present the pictures when k is fixed and let q run from k to infinity. In the left plot of 
Fig. [4j we have k = 1, and we see that the minimal value of \T>\ is above 0.3 and hence 
there is no Fermi surface, even in the approximate sense. In the right plot of Fig. |4j we let 
k = 9, and we find that multiple but finite number of Fermi surfaces emerge. 

We now consider the second major case with r_ < 0. As a concrete example, we choose 
r + = 1 and r_ = —1, and hence we have v = \\Jk 2 + q 2 . Thus the parameters k and q can 
now run independently. We find that for small k or small q, there can be no Fermi surface. 
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Figure 4: These are plots of [D\ with respect to q, for some fixed k. The left plot has k = 1 
and there is no Fermi surface. In the right plot, we have k = 9. Multiple Fermi surfaces 
arise. The black hole parameters are r+ = 2 and r_ = 1/2. 

In Fig. [5j we give such an example with q = 4 and the plot clearly shows the absence of 
a minimum that is close to zero. When k and q are both sufficiently large, Fermi surfaces 
emerge, as can be seen in Fig. [5j 




Figure 5: These are two plots of \T>\ with respect to k, for some fixed q. The left one has 
q = 4 and the values of \T>\ lie roughly in the range of (0.7, 2). The right one has q = 20, 
and multiple Fermi surfaces indeed emerge. The black hole parameters are r+ = 1 = — r_. 

So far we have studied the Green's function with vanishing cj, which allows us to de- 
termine the Fermi surfaces. It is also of great interest to study how the Green's function 
behaves when we change lv. Our analytic expression for general (uj, k) makes the task easy. 
In Fig. [6j we present a graph on how the Green's function behaves under uj on a given Fermi 
surface k = kp. We see that multiple but finite maxima arise for \G(uj, kp)\. 

4.5 Near-extremal limit 

As we have discussed in the previous subsections, there is only one Fermi surface in the 
extreme limit, corresponding to v = 0. For a non-extremal black hole, multiple Fermi 



24 




Figure 6: This is the plot of \G(oj, kp)\ with respect to uj on a Fermi surface. The black hole 
is specified by r+ = 2 and r_ = 1/2. The spinor charge is q = 50. The Fermi surface is at 
kp = 17.998011219, for which the Green's function has a maximum around 10 9 . Multiple 
and finite number of maxima arise as we change u. The maximum in the left is of 10 10 , 
the maximum in the middle with lo = is of order 10 9 , whilst the maximum at the right is 
much smaller of 500. 

surfaces emerge for sufficiently large q. It is of interest to study the near-extremal region. 
Let us redefine the variables r + = r$ and r_ = tq — 5. The near extremal solution can be 
defined as having non-zero 5 but with 5/tq « 1. 

In the literatures, there have been discussions on obtaining decoupling limit of the near- 



extremal black holes. For the black hole (15) with (18), this can be achieved first by the 
following scalings and redefinitions: 

t 
e 



r+ = r + ep , 



r - ep- 



ro + ep. 



t 



(74) 



Then sending the parameter e — > yields the following charged AdS2 x T 2 : 
dp 2 



ds z 



Pi 



{p 2 - p^dt 2 + r£(dx 2 + dy 2 



A = (p- p )dt. 



(75) 



The first two terms appear to describe a two dimensional non-extremal black hole; however, 
it is the AdS2 metric. In order to achieve this AdS2 "black hole", the decoupling limit 
e — >• has to be taken so that the topology R 2 x T 2 of the non-extremal black hole becomes 



AdS2 x T 2 . In this limit, we have r + 



2e/?o — >• 0, and hence the solution should 



really be viewed as the extremal limit, rather than the near-extremal solution. In fact, in 
this topology-changing decoupling limit, all the thermodynamical quantities are identical to 



those of the extremal solution. Thus the solution ( 75 ) should be really viewed as a special 
slice of the AdS2 x T 2 in the extremal limit. In other words, as long as the solution is not 
extremal, no matter how near, the decoupling limit described above does not exist and the 
topology of the black hole remains R 2 x T 2 . 
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Now let us examine how the Green's function (47) behaves in the near-extremal region 
for the oj = case. Let us consider a concrete example, with r + = 1, <5 = 1/10 4 and q = 1. 



For v to be real, we must have k > ko = 3\/ITTT/100. We find that the absolute value |G(A;)| 
lies in between 16.0425 and 1 as k runs from ko to oo, suggesting that there is no Fermi 
surface for real u, as we see in the left figure in Fig. [7j On the other hand, If we let k runs 
to be less than ko such that v becomes pure imaginary, Fermi surfaces emerge, as we can 
see from the right figure in Fig. [7| Thus we see that the characteristics of the Fermi surface 
structure of near-extremal black holes is more or less the same as the generic non-extremal 
black holes. This is consistent with the fact that the topology of the near-extremal solution 
is still a R 2 x T 2 , rather than the extremal AdS2 x T 2 . It is worth pointing out that when 
5 is smaller so that the temperature is lower, the allowed the value of q for the existence of 
a Fermi surface becomes smaller, as one would have expected. 




Figure 7: In both figures, we have u = 0, r+ = 1, r_ = 1 — 10~ 4 and q = 1. In the left 
figure, we plot |Gi| vs. k, with k > kg. In this case v is real and consequently there is 
no Fermi surface. In the right figure, we plot \D\ vs. k, with < k < ko- We have pure 
imaginary v and Fermi surfaces emerge. 



Let us examine the Taylor expansion of the Green's function in terms of 5. Let r + = 
and r_ = ro(l — 5) with 5 « 1, and we keep the linear-order term in 5. Let us first consider 
v is real. 

Case 1: v > \ and v ^ half integers: 



G= 17^ + ^ 5 f S + 0(5 2 ). 

\k-qr 2(A; - qr )i^k + qr {4k 2 - (1 + 4g 2 )r 2 J 

It is advantageous to express the result in terms of v and q variables, given by 
G 



q ± ^/u ^ +q 2 + qr (is 2 + q 2 ± gyV + q 2 ) ^ + 



v 2 {\ - 2i/)r 
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(76) 



(77) 



Case 2: < v < \: 



G 



k + qr 



+ 



k 



_ 5 2u . 



k-qr k- qr 

2~ 4l/ °r(-2z/o)r(l + 2v - 2iq)T(2u + 2ig)(sin(27rg) - i sin(2z/ 7r)) 



1 + 



2vrr(2^ ) 



, (78) 



where i/q = \fk? — y < / 

If instead, we consider v to be imaginary. For a given g, no matter how small, as 5 
becomes smaller, there will be always a Fermi surface. In fact there is an oscillatory factor 
8 U ° in the Green's function, implies that as 5 becomes smaller and smaller, more and more 
Fermi surfaces emerge for k lying between and qr$. The maximum value of k for the Fermi 
surfaces becomes k = qr^ in the extremal limit. As S approaches zero, the number of Fermi 
surfaces for k < qro becomes infinite, signalling instability. This is the physical origin of 
the earlier observation that the wave function near the horizon of the extremal black hole 



for imaginary v has an infinite number of oscillations, as we saw in (71 ). In Fig. |8j we use 
an explicit example to plot this behaviour. 




Figure 8: Here we plot verse log(l/<5). It is clear that |Gi| alternates between and oo 
periodically in terms of log(l/5) as 5 approaches zero. The parameters of the black hole are 
chosen to be r+ = 1, r_ = 1 — 8. The Dirac fermion is chosen to have q = 1 and k = 9/10 
so that v is pure imaginary as 5 becomes small. 



5 Fermi surfaces from general extremal black holes 

As we shall show in section 6, the wave solution for the most general black hole is a Hern 
function whose properties are not conducive for discussing the Green's functions and Fermi 



surfaces in fine detail. In this section we discuss the general extremal black hole (28). The 
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wave functions for general (uj, k) are still of complicated confluent Hern's functions; however, 
they reduce to hypergeometric functions when uj = 0. This is enough to study the existence 
of Fermi surfaces. We find that the general solution of u_ is given by 



ci(l-^)\l- r j : ) l \F 1 [v + i q ,\+v + i q -l + 2v,- ^y ri 
+C2 ( 1 " r {]'\ 1 ~ r iT^- V + *,l->' + W,l- ^ ~ r ^$] • (79) 



where 



k 2 



(80) 



y r (r - n) 

It is clear that for positive v, regularity of the wave function on the horizon r = tq requires 
that C2 = 0. As we discussed in section 4, the ability to fix the boundary condition for u_ 
on the horizon with uj = is the property enjoyed only by the extremal solutions. If v is 
imaginary, the wave function becomes oscillatory when r — > r$ and we shall not consider 
such a case for the reasons that we have discussed in section 4. 

Having obtained u_ , we can easily obtain u + and hence the Green's function, following 
the procedure spelled out in section 3. We have 

Gl = -G 2 ^ = ii^ , (81) 

where 

(y + iq)r 2 F X [\ + v + iq, 1 + v + iq; 1 + 2v; -^] 

7 = 1 — — • (82 

k 2 F 1 \v + iq, \ + v + iq; 1 + 2i/, - ^] 

This expression follows the same form of the Green's function in the earlier examples. In 

fact, the Green's function with uj = for the extremal solution is always real, provided that v 

is real. This is because in the extremal black hole, the solutions for wave functions with real 

v have two branches: one diverges on the horizon and the other converges. The boundary 

condition selects the converging branch for both u+ and n_, which has a consequence that 

they are complex conjugate to each other. For a generic non-extremal solutions, both 

branches for uj = wave solutions converge, and we have to use the w / solutions to 

separate in-falling or outgoing branches and the in-falling solution should be chosen. This 

has a consequence, as we have seen in the non-extremal example discussed in the previous 

section, that u± are not complex conjugate to each other and the Green's function G(Q, k) 

is hence not real. 



The manifestly real expression for the Green's function (79) is given by 

Gl = 1 = ' (83) 
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where 



Z = (l 



1 v 



+ iq) — iff, 5 + v — iq; 1 + 2z^; 



(84) 



It is easy to verify that when q = 0, the Green's function becomes the identity, and hence 
there is no Fermi surface for q = 0. When ri = 0, we recover the Green's function given in 



(63), for which there is only one Fermi surface. 

It is of interest to note that when v = 1/4, the function Z becomes simpler, giving 



Z = (l-i)(l-4ig)(-l-4ig+ a/1 + 16g 2 ) sin ( |(1 + 4ig) arctan 



(85) 

r - ri / 

where we have dropped some inessential real multiplying factors and we have also chosen k = 
r o( r o — r i) for solving v = i. Since z/ = 1/4 corresponds to some Non- Fermi 
liquids, let us analyse the case further. For r\ > 0, there appears to be no Fermi surface. 
However, when T\ < 0, Fermi surfaces emerge. Let us consider a concrete example by 
choosing r\ = —1 and r$ = 16/9. For this choice of parameters, we have k = 5^1 + 16g 2 /9. 
The Green's function is now given by 



(Gi) 2 = G 2 



1-7 



7 



4^(4 - 5 cos(4g log 2)) + 3 sin(4g log 2) 



(86) 



2 4(1 + 7)' 7l6VTl(5-cos(glogl6)) 
Note that 1 — 7 2 is a perfect square and the Green's function is manifestly real. We present 
the plot of the Green's function as a function of q in Fig. [9} It oscillates between 00 and zero, 
corresponding to Fermi surfaces in the standard and alternative quantization respectively. 





Figure 9: For both graphs, we have the extremal black hole with 77 = —1 and tq = 16/9. 



In the left figure, we fix v = 1/4, so that k = 5-y/l + 16<7 2 /9. The plot describes G\ 
with respect to q, and we see it oscillates between infinity and zero, corresponding to 
Fermi surfaces in standard and alternative quantisation respectively. In the right plot, we 
consider the extreme non- Fermi liquids with v = 0, corresponding to k = 20q/9. Fermi 
surfaces emerge also. 

Thus we see that for some appropriately chosen q as indicated in Fig. [9j there exists a 
Fermi surface with v = 1/4. For either the standard or alternative quantization, q must 
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satisfy 1 — = 0, implying 

I2q sin(4g log 2) + 5 cos(4g log 2) - 4 = . (87) 

Unfortunately, there appears to be no analytical solution to this equation, albeit it is rather 
a simple equation. The extreme case of the Non-Fermi liquids corresponds to v = 0. For 
the same black hole, this corresponds to k = 20g/9. We can see in Fig.|9]that such extreme 
situation of Non-Fermi liquids can also arise in this extremal black hole. 

In the above we discussed a special case with v = \, corresponding to some non- Fermi 
liquids. We also saw that the extreme non-Fermi liquids with u = can also arise. By 
numerating examples, we find that it appears that Fermi surfaces with all possible v can 
arise in the (k, q) plane. For the black hole with r\ = — 1 and ro = 16/9, we present the 3D 
plot in Fig. [l0| Comparing to the similar plot in the non-extremal case, we find that the 
Green's function in the extremal background is much simpler than that in the non-extremal 
background. 



Figure 10: The black hole parameters are r\ = —1 and ro = 16/9. This is the 3D plot 
of G in the (k, g)-plane, viewed from the top. Interestingly the Fermi surfaces seem to lie 
naturally in the linear line of k of q, with the diagonal line corresponding to v = 0. 

The situation for the general case is similar. We find no evidence for Fermi surfaces 
when r\ is positive. However, when r\ is none positive, Fermi surfaces emerge. When 
r\ = 0, it reduces to the extremal solution discussed in section 4. In this case, there is 
only one Fermi surface, with k = qr^ for the standard quantization and k = —qro for the 
alternative quantization. When (— r\)/ro « 1, there is still only one Fermi surface in each 
quantization, with k ~ ±q(ro — \r\). The more negative r\ is, the more number of Fermi 
surfaces we can have for each given q. 

For a given black hole with ro and a negative r\ , the existence of Fermi surfaces favors 
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a larger q. We find that there is always at least one Fermi surface and the number of 
the Fermi surfaces increases when q becomes larger. We demonstrate this with a concrete 
example. Let us consider ro = 1 and ro = —1. The Fig. 11 shows two plots of Re(Z) with 
respect to k which runs from \piq to oo. It is clear that for smaller q = 1 in the left figure, 
there is only one Fermi surface and for larger q = 10, there are six Fermi surfaces arising. 




Figure 11: In both graphs, we have ro = 1 and T\ = —1. The plots are Re(Z) with respect to 
k whose minimum value is y/2q for real v. The zeros of the plots indicate the Fermi surfaces. 
The left figure corresponds to q = 1, which only gives rise to a single Fermi surface. The 
right figure corresponds to q = 10, and six Fermi surfaces emerge. For this extremal black 
hole, we find that the minimum Fermi surface has v < ^. 

As we mentioned earlier, when — ti/tq — > + , there is only one Fermi surface with v — > 0. 
If we increase the value of — ri/ro, more Fermi surface will arise. For fixed q, the minimum 
v increases and can eventually pass v = ^. However, for such black holes, we can always 
lower the q value so that the minimal v is less than ^. 



6 Green's function from the most general black holes 



In this section, we study the Dirac equation (29) on the most general black hole ([5]) with 
T 2 topology (e = 0). It is advantageous to use the three roots r%, T2 and = tq of / as the 
parameters of the solution. In this parametrisation, the solution is given by 



/ = i(r-ra)(r-r 2 )(r-ro), A = Q(^- 1 -)dt, 
Q 2 = (nr 2 ) 2 + (Jon) 2 + (r r 2 ) 2 - r 1 r 2 r Q {r 1 +r 2 + r ) 



(88) 



Note that without loss of generality, we have chosen ro > to represent the horizon radius. 
If both r\ and r 2 are complex numbers, they must be conjugate to each other; if they are 
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real, neither of them is bigger than ro- Following the procedure in section 3, we find that 
(u> + qQ/r )r - qQ . „ „ ik 



Ai(r) 



A 2 (r) 



(89) 



(r - n)(r - r 2 )(r - r ) ' rv ^~7^~7^rr7^) ' 

Note that qQ appears in the Dirac equation together, implying that the complicated ex- 



pression for Q in (88) can be absorbed in q such that qQ remains a simple expression. Let 



us first solve for it- . In order to simplify the equation, we let 

( r _ ri )-i(Q+fii)( r _ r2 y(Q-n 2 )( r _ rQ )-Xk> v ( z ) ( 



U- 



(90) 



where 



^2 



n(r - r ) 

r(r - n) ' 
wr 2 



Q 



qQ 



ro{n - r 2 ) ' 



fii 

wr 



(n - r )(ri - r 2 ) ' 



(91) 



(»*2 - r )(r 2 - n) ' "" u (r - n)(r - r 2 ) 4ttT 

Here T is the temperature of the black hole, whose thermodynamics was analysed in section 
2. Note that we have the identity Q\ + f2 2 + £Iq = 0. Thus the three factors (r — rj) with 



{i = 0, 1, 2) in (90) can be simultaneously replaced by (1 — r^/r). Substituting ( |90| ) into the 
wave equation of u_, we find that t>(z) satisfies the following Heun's equation 



v"(z) + 



|-2igo |-2i(Q+Oi) |+2i(Q-n 2 ) 
z 2—1 2— a 



z(z — l)(z — a) 



v(z) = 0. 



with 



(r - r 2 )n 



A' 2 



(92) 



(93) 



(r - ri)r 2 ' (r - ri)r 2 ' 

Therefore, the wave function u_ can be solved in terms of the following general Heun's 
function Hi: 

u _ = Cl ( r _ ri )-KQ+ni)( r _ r . 2 )i(Q-fi2)( r _ r . o )-in 0x 
i»?(a, &; 0, |, \ - 2ift , g - 2i(Q + Qi); -z) 
+c 2 r-i~ 2iQ °(r - n)~ i(| 3 +Ql) (r - r 2 ) i( ^~ fi2) (r - r )^ iS! °x 
flgfo, /3; | + 2ifi , 1 + 2iO , 1 + 2iO , \ - 2i(Q + fii); , (94) 



where 



/3 



+ 2i^ ) 



2igQ 



+ 2ri : 2 : (ri + r2)ro (i + 4i^ ; 



(95) 



' """ / Vr 2 (ro - ri) ' 2r 2 (r -ri) 
As we have mentioned in section 4.4, the parameter combination u, which plays an impor- 
tant role in extremal solutions, has no apparent physical significance in the non-extremal 
backgrounds. In fact it does not appear in the above wave solutions at all. Comparing to 
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(80) for the general extremal case corresponding to having r<i = ro, we expect that the v 



may arise in the following combination 



Q 2 ) 



T2— >ro 



(96) 



—¥ extremal 

Thus we see clearly that in the non-extremal background, such a combination does not 
appear in the solution and hence it should play no particular role. There is no physical 
reason to impose the reality condition. It follows that there is no constraint on (k, q) in the 
non-extremal case. The emergence of the z^-like quantity in the subclass of non-extremal 
black holes discussed in section 4 was pure accidental. 



On the horizon r = ro, corresponding to z = 0, the Heun's functions in (94) all reduce 



to identity, and hence the c\ term corresponds to the in-falling mode and the ci term is 
the outgoing mode. For the black hole background, it is necessary to choose the in-falling 
solution only, and hence we set c<i = 0. Having obtained the in-falling solution for u_, the 



function u + can then be read off from (39). (It actually requires some derivative identity 
given presently.) Thus the full set of in-falling solutions are 

u _ — (l _ nj-i(Q+f^i) ^ _ r2_y(Q-n 2 ) ^ _ r^^-iQo x 



H£[ a, b; 0, A, ± - 2ift , 5 - 2i(Q + Qi); -z 



u + 



■in . 



Hl[a,F;\ -2ifi ,l 



2i^o> § 



2ifi ,| +2i(Q-Oi); 



where (5* is the complex conjugate of (3 given in (95) and c is a complex constant: 



k 



ro > 



ro - 



(97) 



(98) 



This constant coefficient can only be obtained from evaluating (39), even though the struc- 



ture of the solution for u+ can be obtained right away by the fact that u+ takes the complex 



conjugate form of u_. It follows from (40) that the Green's function can now be obtained 
straightforwardly: 



G(u,k) = i 



1 --y(u,k) 



where 



j(u,k) 



l + 7(w,fc) ' 

cHi(a,f3*; \ - 2ifi , 1 - 2ift , \ - 2i^o, \ + 2i(Q - fii); 



r -n 



H£{ a, b; 0, A, \ - 2ift , | - 2i(Q + 



(99) 



(100) 



Compared to the hypergeometric functions, Heun's functions are much less studied. In 
particular, in the Mathematica package, the Heun's functions are not yet coded. In the 
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Maple package, although the function is built in, very limited properties are coded. Note 
that when r\ or r 2 vanishes, the Heun's functions are reduced to hypergeometric functions 
and we obtain the results for the "BPS" black holes, discussed in section 4. If we let t\ or 
r 2 to be ro, we obtain the results for the most general extremal solution. The limit is rather 
subtle since several arguments in the Heun's function diverge. We find that the Heun's 
functions become some confluent Heun's functions Hie in this limit. The proper procedure 



requires the identity ( 104 ) with precise coefficients which are yet unknown in the literature. 



We present the general solution for u_ for general extremal solution for which r 2 = r$: 
n_ = c 2 (l - a)-|+i(Q-n) (1 _ myi(Q-a) H£c ^ _ _i ; _x_ + 2i( g _ ^ ^ _^ 

+ci(l - r -xy(Q-n) {1 _ a)-i(Q-6) Mc ( _ 2 m,-l -\ + 2i(Q -Q),5, 77; -z) ,(101) 

where 

q _ qQ ^ _ urx ^ _ (rp - n)r 

roiro-ri)' (r - n) 2 ' ri(r-r )' 

^ ^ A; 2 
S = iO(l + 2i(Q - U) , r/=|-iiQ , -, (102) 

In order to proceed further, it is necessary to fix the horizon condition on this general 
wave solution for the two-parameter extremal black hole. However, on the horizon, we have 
z — > 00, it is thus necessary to change the variable z to 1/z. Unfortunately, we have not 
found in the literature the necessary identities to perform this transformation. We leave 
this problem open for future investigation. 



In order to determine the Fermi surfaces from (99), we set oj = 0, and the Green's 
function becomes 



G(0,k) 



l + 7(*0 ' 

i(r - n)i- 2i « (r - r 2 )§+ 2i( ? Hi(a, f, 1, |, \ + 2iQ; 
l(k) = 7 — 4 = \ ^,(103) 



kr He(a,b;0±±±-2iQ; 



, 2 , 2 , 2 ^xvg, rQ _ ri 



where /3 is given in (95) with u; = 0. Analogous to the simpler non-extremal example 
discussed in section 4, G(0, fc) is not real and there is no literal divergence . However, there 
are effective divergencies corresponding large spiked maxima of |Cr(0, k)\. To demonstrate 
that the above Green's function can indeed produce effective Fermi surfaces, let us consider 
a concrete example. Let us consider ro = 3, T\ = 1 and r 2 = —2. A fermi surface occurs 
at qQ = 11.75535749 and k F = 0.7903492261. The absolute value of the Green's function 



at this point is of order 10 9 . We present some graphs of |G(0, k)\ in Fig. 12. It is easy 
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to see that the maxima associated with the effective Fermi surface is "spiked", in that for 
8k ~ 1 around the kp, we have S\G\ ~ 10 9 . We also consider the black hole with ro = 3, 
r\ = 1 and ri = 2. In other words, all the roots are positive. We find a Fermi surface at 
qQ = 5.3042024 and kp = 0.264110954. In this case, the spiked maximum of G\ is also in 
the order of 10 9 . 





Figure 12: For both plots, the black hole parameters are ro = 3, n = 1 and ri = —2. The 
left is the plot of \G\\ with respect to k, with qQ = 11.75535749. The maximum occurs at 
k F = 0.7903492261 and it is of order 10 9 . The minimum is of order 10 8 . The right is the 
three-dimensional plot of \G\\ with both k and qQ variables. Many spiked maxima arise. 



The final example is provided in Fig. 13 on the black hole ro = 3, r\ = 1+i and r2 = 1— i. 
In this case, / has only one real root, and hence there is no inner horizon. Again, the 3D 
figure of |G| in the (k, qQ)-plane shows a fascinating rich structure of spiked maxima. 




Figure 13: The 3D plot of |G| in the (k, qQ)-p\ane shows a fascinating rich structure of 
spiked maxima. The black hole is specified by ro = 3, r± = 1 + i and T2 = 1 — i. 



To end this section, We are obliged to acknowledge that we have used extensively the 
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properties of the general Heun's functions listed in |3 1 1. 132j. For our purpose, we find that 
the following identities are particular useful. 

(1) The Heun's function H£(a,b; a, P,j,5; z) can be written as the linear combination 
of the following two Heun's functions: 

z- a M(i,a(e-/3) + f(5-/3) + ^;a,a- 7 + l,a-/3 + l,5;i), 
z-f>Bt(i,P{e - P) + f (5 - P) + b -- p, - 7 + 1, P ~ a + 1, 5; i) . (104) 

Unfortunately, the precise coefficients are not yet known in the literature. Owing to the 
lacking of this knowledge, we were unable to perform the subtle extremal limit of the wave 
solutions and the Green's function from those of most general non-extremal backgrounds. 

(2) The following is a derivative identity that we find useful to calculate u + from u_ : 

^H£(a, b; 0, p, 7, 5; z) = H£(a, b + 7 + (a - 1)5 + P + 1; 2, P + 1, 7 + 1, 5 + 1, z). (105) 

(3) We also used the following identity to put itl? in the form as the complex conjugate 
of u° ut : 

H£(a, b; a, p, 7, 5; z) = (1 - ^(l - |)"^+7+<5 x 

m(a, b - y((5 - l)a + p - 7 - 5); 7 - P + 1, 7 - a + 1, 7, 2 - 5; z) .(106) 

7 Conclusions 

In this paper, we considered the conformal gravity in four dimensions. It is constructed from 
the Weyl-squared term together with the minimally-coupled Maxwell field. The most gen- 
eral spherically-symmetric charged AdS black holes up to conformal factor was constructed 
in |15j . We generalize the result to the torus T 2 and hyperbolic H 2 topologies. We analyze 
the global structure and obtain the full set of thermodynamical quantities that satisfy the 
first law of thermodynamics. 

Aside from the facts that conformal gravity plays an important role in critical gravity 
and also in M = 1, D = 4 off-shell supergravity, the main reason that we construct the 
T 2 -symmetric charged AdS black holes is to use them as gravitational backgrounds to study 
the properties of the boundary field theory. Charged AdS black holes such as the extremal 
Reissner-Nordstrom solutions in general dimensions were demonstrated to be dual to Fermi 
and non-Fermi liquids and Fermi surfaces were determined. In the literature, the analytic 
Green's functions for general (w, k) are lacking and properties of G(u>, k) were obtained 
through numerical analysis. In this paper, we have demonstrated that the Dirac equation 
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for a charged massless spinor on the AdS black holes in conformal gravity can be solved 
exactly. This allows us to obtain the exact Green's function for all (oj, k). Since in conformal 
gravity the black holes are given up to arbitrary conformal transformations, this leads to 
the question how general our results are. We show that the Green's function, and hence 
the Fermi surfaces are invariant under the conformal transformations that preserve the T 2 
isometry. Thus our results cover all the T 2 -symmetric black holes. Some salient results of 
this paper were reported earlier in |30| . 

The Green's function from the general black hole is expressed in terms of the general 
Heun's functions. The Heun's functions are much less studied and there is no package in 
Mathematica and only a limited package in Maple. This makes it difficult to analyse the 
Fermi surfaces and their properties in detail using these programs. Nevertheless, we made 
numerical plots of (7(0, A;), which show fascinating rich structures of spiked maxima in the 
plane spanned by (k,q). 

The most general black hole solution contains three non-trivial parameters. We consider 
some special subclasses of the black holes with two non-trivial parameters: the horizon 
radius and the charge. The solution can be viewed as being pseudo-supersymmetric from 
the point of view j\f = 1, D = 4 off-shell supergravity and was obtained in [16]. We find 
that the Heun's functions are reduced to hypergeometric functions. The Green's function 
for general (u>, k) becomes fully analysable. The Fermi surfaces are defined by the poles of 
G(0, k). In the extremal limit of this subclass of backgrounds, we find that there can only 
be one Fermi surface, corresponding to having v = 0. We obtain the small uj behaviour 
near the Fermi surface. When the background is non-extremal, there is no Fermi surface at 
all for real v. However, Fermi surfaces emerge when v is pure imaginary. 

Our explicit results demonstrate that there are qualitative differences between the ex- 
tremal background whose decoupling limit is AdS2 x T 2 and the non-extremal background 
whose topology is R 2 x T 2 . For both extremal and non-extremal backgrounds, if we have 
wave solutions for generic (uj, k), we can organise the general solution in terms of in- falling 
and outgoing modes on the horizon. The boundary condition of a black hole horizon selects 
the in-falling solution, which is associated with the retarded Green's function. Since the 
Fermi surfaces are determined by the Green's function with ui = 0. It is thus tempting to 
think that we only need solve for wave functions with uj = 0, in which case, the equations 
tend to be reduced significantly. However, once u is set to zero, we cannot organise the 
solutions in terms of the in-falling and outgoing modes. In the extremal background, for 
real u, the boundary condition on the horizon could still rule out one of the two modes, since 
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it diverges on the horizon. Furthermore, if v is imaginary, the modes become oscillatory in 
the near-horizon region. This signals instability of the system. 

The situation is quite different in the non-extremal case. When ui is set to zero, we find 
that both modes are finite on the horizon. Thus for non-extremal backgrounds, there is no 
procedure to select the correct mode if we have only solutions with co = 0. This implies 
that in order to obtain the Fermi surfaces, we have to obtain solutions with non-vanishing 
lj, even though the Fermi surfaces are determined by the properties of the Green's function 
with vanishing cu. Another difference is that the in- falling wave solution with uj = is not 
oscillatory even when v is imaginary. Thus there is no instability associated with imaginary 
v. Using our explicit example, we show that for imaginary v, the number of Fermi surfaces 
increases as the black hole approaches the extremality, and the number becomes infinite in 
the extremal limit. This provides an explanation of the origin of the instability associated 
with imaginary u in the extremal backgrounds. 

We also study the Green's function for the general extremal solution with two non-trivial 
parameters, characterised also by the horizon radius and the charge. The Green's function 
G(oj, k) is now expressed in terms of confluent Heun's functions. However, G(0, k) is given 
in some simpler hypergeometric functions. We find multiple numbers of Fermi surfaces 
arising in this case. These Fermi surfaces can be ones for either Fermi liquids (y > 1/2) or 
non- Fermi liquids (y < 1/2), including the extreme case v = 0. 

Our further analysis on the wave solutions and the Green's function for the most general 
black hole indicate the v parameter that appears in the non-extremal "BPS" solution is 
purely accidental. In the Green's function and the wave solutions in the most general black 
hole, there is no parameter combination that resembles v. The reality condition on v for 
the extremal background which implies k > k m i n (q) is relaxed in the non-extremal case 
where (k, q) can take all real values. The numerical plots of the Green's function in the 
(k, q) plane show a fascinating rich structure of spiked maxima. 

Conformal gravity that involves the Maxwell field is uniquely four dimensional and the 
dual field theories have (2 + 1) dimensions. Our results show that conformal gravity and its 
charged black holes provide useful playgrounds for investigating properties in both Fermi 
and non-Fermi liquids on a two-dimensional plane, not only at zero temperature but also at 
any finite temperature. The fact that the Dirac equation for a charged massless fermion is 
invariant under the T 2 -preserving conformal transformations implies that conformal gravity 
is particularly suitable for studying such fermionic system. 
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